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We review a recently proven Lieb-Robinson bound for general, many-body quantum 
systems with bounded interactions. Several basic examples are discussed as well as the 



■ 1. Introduction 

o 



Much physical intuition is based on locality properties of the system under consid- 
. eration. Objects (or particles) are associated with regions (or points) in space and 

non-trivial interactions typically occur over short distances. For systems governed 
by a relativistic dynamics, the time evolution preserves this notion of locality. 

In condensed matter physics, however, many interesting physically relevant phe- 
nomena are modeled by quantum many-body systems, e.g., super-conductivity and 
Bose-Einstein condensation. More abstractly, the models of quantum computation 
and discrete versions of field theory are described in terms of quantum lattice sys- 
' terns. For these non-relativistic systems, defined e.g. by a Hamiltonian with nearest 

neighbor interactions, the associated dynamics does not preserve locality in the 
sense that there is no strict equivalent to a finite speed of light. 

Remarkably, Lieb and Robinson, see [1], proved that an approximate form of 
locality, which we refer to as quasi-locality, does hold for the dynamics associated 
to certain quantum spin systems. This important result establishes the existence 
and proves a bound for an approximate light cone which limits the rate at which 
disturbances, as evolved by the dynamics, can propagate through the system. More 
concretely, they proved that a local observable evolved for a time t > remains 
essentially localized to a region of space whose diameter is proportional to t. They 
dubbed their estimate as a bound on the system's group velocity, but we prefer 
to describe the analysis as a Lieb-Robinson bound and the resulting estimate: the 
Lieb-Robinson velocity. 

After the initial result by Lieb and Robinson in 1972 and some calculations for 
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specific models [2] a few years later, these locality estimates for quantum systems 
received relatively little attention. It was not until Hastings' impressive work of 
2004, see [3], that a genuinely renewed interest in these bounds was established. 
Since then, a number of generalizations of the original result [4-11] and a wealth 
of interesting applications [3, 12-27] have demonstrated the importance of these 
bounds. 

In this brief note, we introduce the general set-up for Lieb-Robinson bounds in 
the context of quantum systems with bounded interactions. We describe the cor- 
respondence between quasi-locality and the usual commutator estimates typically 
referred to as Lieb-Robinson bounds. The final section is a short list of further 
generalizations and a variety of applications. 

Before we begin, we make the following useful observation. In quantum mechan- 
ics, one is often interested in a single quantum system, i.e., a specific Hilbert space 
H and a densely defined, self-adjoint operator H . For each normalized vector, or 
state ip £ H, the solution of the Schrodinger equation 



governs the dynamics of this system. The solution, of course, is t/j(t) — e~ ltH ij}. 
For particles in a domain A C W , the typical Hilbert space is H = L 2 (A) and it 
is common to have H a self-adjoint realization of the Laplacian. Corresponding to 
any normalized vector ip £ L 2 (A) the solution tp(t) is called a wave function, and it 
is interesting to investigate the evolution of this function in space, i.e. A. 

This is not the form of locality established by a Lieb-Robinson bound. Lieb- 
Robinson bounds are about collections of interacting quantum systems distributed 
in space. The bounds estimate the rate at which disturbances propagate through 
this collection. 

2. The General Set-Up 

As discussed above, Lieb-Robinson bounds estimate the rate at which disturbances 
propagate through a collection of quantum systems. The basic set-up is as follows. 

2.1. Collections of Quantum Systems 

Let T be a countable set, and consider a collection of quantum systems labeled by 
x £ r. By this, we mean that corresponding to each site x £ T there is a Hilbert 
space T-L x and a densely defined, self-adjoint operator H x acting on H x . The operator 
H x is typically referred to as the on-site Hamiltonian. For finite A C T, the Hilbert 
space of states corresponding to A is given by 



(1) 



H A = (g)H : 



X ) 



(2) 



and the algebra of observables is 



Aa = ®B(H x ) = B(H a ) 



(3) 



xeA 



November 23, 2010 1:30 



WSPC - Proceedings Trim Size: 9.75in x 6.5in q'mathllsims 



3 

where 13(H) is the set of bounded linear operators over the Hilbert space H. Thus 
an observable A e A a depends only on those degrees of freedom in A. Of course, 
for any hnitc X C A, an observable A e Ax can be uniquely identified with the 
observable A <g> 1a\x £ -4a > an d therefore Ax C Aa- 

In general, these collections of quantum systems are used to describe many inter- 
esting physical phenomena e.g., the moments associated with atoms in a magnetic 
material, a lattice of coupled oscillators, or an array of qubits in which quantum 
information is stored. Below we indicate two important types of examples. 

Example 2.1. A quantum spin system over T is defined by associating a finite 
dimensional Hilbert space to each site x € T, e.g., H x = C" x for some integer 
n x > 2. The dimension of H x is related to the spin at site x by n x = 2J X + 1, 
i.e. n x — 2 corresponds to spin J x = 1/2, n x = 3 corresponds to spin J x = 1, 
etc. As an on-site Hamiltonian, a common choice is to select a spin matrix in the 
n^-dimensional irreducible representation of su(2). When n x = 2, these are just the 
Pauli spin matrices: 

Example 2.2. A quantum oscillator system over T corresponds to associating an 
L 2 space to each site of V, e.g. one can take H x = L 2 (M.) for each x G T. In contrast 
to the previous example, each single site Hilbert space is infinite dimensional, and 
moreover, the on-site Hamiltonians are typically functions of position q x , the multi- 
plication operator by q x in L 2 (R, dq x ), and momentum p x = — i-^-; both unbounded 
self-adjoint operators. 

Despite the fact that these examples are quite different, the general techniques 
described below apply equally well to both cases. 

2.2. Interactions and Models 

The systems described above are of particular interest when they are allowed to 
interact. In general, a bounded interaction for such quantum systems is a mapping 
$ from the set of finite subsets of T into the algebra of observables which satisfies 

&(X)* = <$>(X) e A x for all finite Icf (5) 

A model is defined by the set T, the collection of quantum systems {(H x , H x )} xer , 
and an interaction $. 

Associated to a given model there is a family of local Hamiltonians, {i^A}, 
parametrized by the finite subsets of T. In fact, to each finite A C T, 

H A = £ H x + £ $(A) (6) 

xeA XcA 

is a densely defined, self-adjoint operator. Here the second sum is over all finite sub- 
sets of A, and is therefore finite. By Stone's theorem, the corresponding Heisenberg 
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dynamics, r t A , given by 

t^{A) = e itHA Ae- itHA for all A G A\ and t G R, (7) 

is a well-defined, one-parameter group of automorphisms on ,4a- Here are two typical 
models of interest. 

Example 2.3. Fix an integer v > 1 and let T = . Consider the quantum spin 
system obtained by setting T~L X = C 2 for all x G 7L V . Take as on-site Hamiltonian 
H x = S 3 using the notation from (U above. Let 5> be the interaction defined by 
setting 

(SlSl + S*St> + S*S* HX = {x,y} a nd\x-y\ = l, 
*W = \ otherwise (8) 

where for each zeZ", any k G {1, 2, 3}, and each finite volume AcZ", the quantity 

S k z = 1 <8> • • ■ ® 1 ® 5 fc ® 1 (8 • • • ® 1 (9) 

where £ , again from (j4|), appears in the z-th factor of .4a = £>(C 2 ). 
A nearest neighbor, spin 1/2 Heisenberg model on % v corresponds to 

H A = h H x + J *P0 

i£A XcA 

= + j (S^Sy + S^Sy + S^Sy) (K)) 

X£A x,y£A: 

\x-y\ = l 

for all finite subsets A C Z" ■ Here h and J are real- valued parameters of the model. 

Example 2.4. Fix an integer v > 1 and let r = Z". Consider the quantum oscil- 
lator system obtained by setting H x = L 2 (M.) for all x G V . A nearest neighbor, 
anharmonic model on Z" is defined analogously, e.g. with 

H A = Y / Pl + V(q x )+ ^x-q y ) (11) 

l£A x,y€A: 

\x-y\ = l 

for all finite subsets AcT. The parameters of this model are V and <&. Of course, 
V must be chosen so that the on-site Hamiltonian H x = p x + V(q x ) is self-adjoint, 
and $ is assumed to be in L°°(R). 



2.3. Observables and Support 

The support of an observable plays a crucial role in Lieb-Robinson bounds. We 
introduce this notion here. Let T be a countable set and {(T-i x , ffi)}igr a collection 
of quantum systems. As we have seen above, for any two finite sets A G A G T, 
each A G Aa can be identified with a unique element A ® 1a\a G -4a- For this 
reason, A\ C .4 a for any Ao C A. 

Given an observable A G ,4a, we say that A is supported in X C A if A can be 
written as A = A (g) 1a\x with A G -4x- The support of an observable .A is then 
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the minimal set X such that A is supported in X. We will denote the support of 
an observable A by supp(^4). 

Due to the fact that we are considering non-relativistic systems, i.e., models for 
which there is no strict equivalent to a finite speed of light, the following observation 
is generally true. Let X C A C T. Consider a local Hamiltonian H\ defined in terms 
of a non-trivial interaction; assume e.g. the interaction is nearest neighbor. Then, 
for general A £ Ax, s ^PP( T o(A)) = supp(^4) C X, however, supp(r t A (^4)) = A for 
all t ^ 0. Hence, a strict notion of locality, implicitly defined here in terms of the 
support of an observable, is not generally preserved by the Heiscnberg dynamics. 

Lieb-Robinson bounds address the following simple question: Does the Heiscn- 
berg dynamics corresponding to, e.g. short range interactions, satisfy some weaker 
form of locality? 

2.4. From Locality to Commutators 

Lieb-Robinson bounds are often expressed in terms of commutator estimates. The 
relationship between these estimates and the support of local observables is due 
mainly to the tensor product structure of the observable algebras. We briefly discuss 
this fact in this section. 

Let T be countable set and {{V. x , H x )} x£ r denote a collection of quantum sys- 
tems. Consider two finite sets X, Y C T. If A £ Ax, B £ Ay, and X n Y = 0, then 
for any finite set A C T for which lUYcAwe can regard A,B £ A \ and as such 
[A, B] = due to the structure of the tensor product. In words, observables with 
disjoint supports commute. 

Conversely, Schur's lemma demonstrates the following. If A £ Aa and 

[A, t A \ Y (g>B}=0 for all B £ Ay , (12) 

then supp(^l) C A \ Y. In fact, a more general statement is true. If A £ Aa almost 
commutes with all B £ Ay, then A is approximately supported in A \ Y. The 
following lemma appears in [28]. 

Lemma 2.1. Let Hi and H2 be Hilbert spaces and A £ B{H\ 8%). Suppose there 
exists e > for which 

\\\A,-h®B]\\ < e\\B\\ for all B e B(H 2 ). (13) 

Then, there exists A' £ B(T-Li), such that 

\\A'®1 2 -A\\<e. (14) 

Equipped with this lemma, we see that uniform estimates on commutators 
provide approximate information on the support of observables. Thus, for any 
A £ Ax, we can approximate supp(r t A (A)) by bounding || [t^(A), B] \\ for all B 
with supp(B) C Y. Here the estimates will, of course, depend on the distance be- 
tween X and Y in T and the time t for which the observable has been evolved. This 
is the basic idea of a Lieb-Robinson bound. 
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3. Lieb-Robinson Bounds 

In this section, we describe in detail assumptions on the set T and the interactions 
$ under which one can prove a Lieb-Robinson bound. We also present a precise 
statement of the estimate and discuss several relevant consequences. For a proof in 
this setting, we refer the interested reader to [9]. 

3.1. On the Geometry of T 

For many models, e.g both Example 12.31 and !2. 41 the set V = 7L V . In general, though, 
the lattice structure of Z" is not necessary to prove a Lieb-Robinson bound. The 
following assumptions are sufficient. Let T be a set equipped with a metric d. If T 
has infinite cardinality, we further assume that there is a non-increasing function 
F : [0, oo) — > (0, oo) which satisfies two conditions. First, we assume that F is 
uniformly integrable over T, i.e., 



Next, we assume there exists C > such that the following convolution condition 
is satisfied: for all x, y € T, 



The inequality (|16|) is quite useful in the iteration scheme which is at the heart of 
proving a Lieb-Robinson bound. 

Here is an important observation. Let T be a set with a metric and F satisfy 
the properties mentioned above with respect to T. In this case, the function F a 
defined by setting F a (r) = e~ ar F(r) for any a > also satisfies (fT5|) and (fT6|) 
above with \\F a \\ < \\F\\ and C a < C. The choice of an exponential weight here 
is convenient, but not necessary. In fact, G — wF satisfies (|15p and (fl6|) for any 
positive, non-increasing, logarithmically super-additive weight w, i.e. a function w 
for which w(x + y) > w{x)w(y). 

Example 3.1. Consider the case of T = Z v . For any e > 0, the function F(r) = 
(1 + r)~ ( - I/+£ -' is positive, non-increasing, and 



A short calculation shows that the convolution constant for this F satisfies C < 
2^+c+i||^7 , ||_ Thus, such functions do exist. As a final remark, we note that the 
exponential function does not satisfy the convolution condition (1161) on Z", however, 
F a (r) = e- ar /(l + r) v+e certainly does. 




(15) 




(16) 




(17) 
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3.2. Assumptions on the Interaction $ 

Locality estimates are valid when the interactions are sufficiently short range. For 
general sets T, a sufficient decay assumption can be made precise in terms of the F 
function introduced in the previous sub-section. 

Let r be a set with a metric d and a function F as in Section 137X1 For any a > 0, 
denote by B a (T) the set of all those interactions $ for which 

- ^i&fcm £ mxn * °°- (18) 

x,y£X 

Example 3.2. Many interesting models have finite range interactions, see e.g. Ex- 
ample 12.31 and 12.41 An interaction $ is said to be of finite range if there exists a 
number R > for which = whenever the diameter of X exceeds R. In the 

case of r = and F(r) = (1 + r)~ l/ ~ e , it is easy to see that all uniformly bounded, 
finite range interactions satisfy ||$|| a < oo for all a > 0. 

Example 3.3. Another important class of models involve pair interactions. An 
interaction $ is called a pair interaction if $(A) = unless X = {x,y} for some 
points x, y 6 T. In the case of T = V and F(r) = (1 + r)~ l/ ~ e , it is easy to see that 
all uniformly bounded, pair interactions that decay exponentially in \x — y\ satisfy 
||$ || a < oo for some a > 0. In fact, if the pair interactions decay faster than an 
appropriate inverse polynomial, then ||$||o < oo, and this is sufficient for a decay 
estimate on the relevant commutators. 



3.3. The Main Result 

We can now state a Lieb- Robinson bound, proven in [9], for the systems introduced 
above. 

Theorem 3.1. Let T be a set with a metric d and a function F as described in 
Section \3.1i Fix a collection of quantum systems {(H x , H x )} xe r over T, and for 
any a > 0, let $ G B a (T). Then, the model defined by T, {(Ji x , H x )} xe r, and $ 
satisfies a Lieb-Robinson bound. In fact, for each fixed finite subsets X, Y C T , and 
any finite Ac! with X U Y C A, the estimate 



[r t A (A),B]||<2||A||||i?||min 



l,g a (t) F a{d{x,y)) 



(19) 



holds for any A £ Ax, B £ Ay, and t G M. Here the function g a is given by 



(c^(enn a c M _ lhfxnY = t 

5a(t) ~\ C-V"*!^!*! otherwise. {M > 

As a corollary, a more familiar form of the Lieb-Robinson bound can be expressed 
in terms of d(X, Y) = inin xe x,yeY d(x, y), namely 
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Corollary 3.1. Given the assumptions of Theorem \S.l\ and a > 0, the estimate 



[rHA),B}\\< 



2||A|1[|B||||F 



li min[|X|, |y|] e -«(dpW-"*(«)|f|) (21) 



C a 

is valid. Here \X\ denotes the cardinality of X and the quantity v$(a) is given by 

v$(a) = . (22) 

a 

Let us make a few remarks to help interpret these bounds. As we indicated in 
Section F2.2[ the Heisenberg dynamics, r t A , forms a one-parameter group of automor- 
phisms on Aa and so the estimate || [t^(A), B] \\ < 2||A||||B|| is always true. What 
we see from ([T5)| . and more directly in (|2"Tj) . is that if A and B have disjoint supports 
X and Y respectively, then [to (A) , B] — [A, B] — and the estimate on [r f (A) , B] 
is small in the distance d(X, Y) for times 

l«l<*^. (23) 

For this reason, the quantity u$ (a) is called a Lieb- Robinson velocity for the model 
under consideration. In fact, using Lemma [2.1[ we see that for each A € Ax the time 
evolution t^(A) is approximately supported in a ball of radius w$(a)|i| about X. 
Thus the dynamics of the system remain essentially confined to a light cone defined 
by this Lieb- Robinson velocity. Moreover, the velocity u$(a), see (|22|) . which governs 
the rate at which disturbances propagate through the system, depends only on the 
interaction $ and the geometry of T; specifically, it is independent of the on-site 
Hamiltonians. 

Another crucial fact about these Lieb- Robinson bounds is that the explicit esti- 
mates, in particular the velocity, do not depend on the finite volume A on which the 
dynamics is defined. This suggests, and can be proven in this setting see e.g. [27], 
that a thermodynamic limit for the dynamics exists. It too satisfies the same Lieb- 
Robinson bound. 

As has been useful in a variety of applications, it is interesting to note the 
dependence of these bounds on the support of the corresponding observables. Since 
only the minimum cardinality appears, one of two the observables could be allowed 
to be volume, i.e. A, dependent without sacrificing the bound. In fact, a more 
detailed analysis shows that only the minimum cardinality of the boundary of the 
supports of the observables is relevant, see [9] for a precise statement. 

When $ e <8a(r) for some a > 0, the Lieb- Robinson bounds decay exponentially 
in the distance between the supports of the observables. The rate of this exponential 
decay, here the number a > 0, is usually of little consequence. For this reason, if 
$ G £>a(r) for all a 6 (a, f3), the optimal Lieb- Robinson velocity is given by 

2||$||aCa .... 

inf v^(a) — mf . (24) 

ae(a,P) a£(a,/3) a 
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We now estimate this optimum velocity for the previously mentioned examples. 

Example 3.4. Let T = V , F(r) = (1 + r)~ l/_e , and consider the spin 1/2 Heisen- 
berg model introduced in Example 12.31 Clearly, for any a > 0, 



As we observed above, this estimate on the velocity is independent of the on-site 
parameter h. 

Example 3.5. Let L = Z", F(r) = (1 + r)~ v ~ e , and consider the anharmonic 
Hamiltonian introduced in Example 12.41 A similar calculation gives a bound on the 
optimal velocity of 



which is independent of the on-site function V, so long as the self-adjointness as- 
sumption is satisfied. 

4. Some Words on Generalizations and Applications 

Over the past few years applications have driven a number of interesting gener- 
alizations of the original Lieb-Robinson bounds. Several review articles have been 
devoted to many of these specific applications, see [8, 29], and some lecture notes 
from schools on topics concerning locality are now available [30, 31]. In this short 
note, we make no attempt to give an exhaustive list of generalizations and ap- 
plications, but rather we list many relevant works to give the interested reader a 
reasonable starting point to further investigate this active area of research. 

4.1. On Generalizations 

The Lieb-Robinson bound stated in Theorem 13.11 and proven in [9] , already in- 
cludes several generalizations of the original result. Most importantly, it applies to 
quantum systems with infinite dimensional, single site Hilbert spaces. In addition, 
no assumption on the lattice structure of L is necessary, and the dependence of the 
bound on the support of the observables has been refined. 

Recently, Lieb-Robinson bounds have been proven for time-dependent interac- 
tions, see [28]. Moreover, Poulin demonstrated in [11] that these estimates also hold 
for an irreversible, semi-group dynamics generated by Lindblad operators. 

Quite some time ago, it was proven in [32] that the analogue of Lieb-Robinson 
bounds hold for the non-relativistic dynamics corresponding to classical Hamilto- 
nian systems. In the past few years, further work in this direction has appeared 
in [33] and [34]. 



\\®\\ a = e a 2 v+ *ZJ <oo, 
and therefore, a bound on the optimal velocity of this model is given by 



(25) 




(26) 




(27) 
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An important open question is: To what extent do Lieb-Robinson bounds apply 
in the case of unbounded interactions? For certain simple systems, there has been 
some progress on this issue. Lieb-Robinson bounds for general harmonic systems 
first appeared in [7]. It was proven in [9], see also [21, 27], that these estimates 
also hold for anharmonic systems if the perturbation is sufficiently weak. A recent 
result in [35] suggests that such bounds apply much more generally. Finally, a Lieb- 
Robinson estimate for commutator bounded operators appeared in [10]. 

4.2. On Applications 

Many of the generalizations mentioned above came about by pursuing concrete 
applications. As discussed in the main text, the resurgence of interest in Lieb- 
Robinson bounds was mainly motivated by Hastings' 2004 paper [3] on a proof of the 
multi-dimensional Lieb-Schultz-Mattis theorem. In this incredibly influential paper, 
Hastings discussed generalized Lieb-Robinson bounds, an Exponential Clustering 
theorem, and pioneered his notion of a quasi-adiabatic evolution. This single work 
inspired a flurry of activity which continues to this day. 

The Lieb-Schultz-Mattis theorem, see [36], concerns the spectral gap between 
the ground state energy and that of the first excited state for the nearest-neighbor, 
spin 1/2 Heisenberg model in one dimension. They proved that for a finite volume 
of size L, if the ground state is unique, then the gap is bounded by C/L, for some 
constant C. Further generalizations, to models with arbitrary half-integer spin and 
to a statement valid in the thermodynamic limit appeared in [37]. Hastings paper [3] 
developed a multi-dimensional analogue of this result. In fact, his argument yields a 
gap estimate applicable in a great deal of generality, see [15] for a precise statement. 
Recent reviews of these results appear in [8] and [30]. 

The Exponential Clustering theorem is a proof that the ground state expec- 
tations of gapped systems decay exponentially in space. Proofs of this result first 
appeared in [4] and [5]. A refinement of the dependence of the estimates on the 
support of the observables was proven in [8] , and this fact was later used by Matsui 
in [26] to investigate a split property for quantum spin chains. 

It is well known that, for quantum spin systems, a Lieb-Robinson bound may be 
used to establish the existence of a thermodynamic limit for the Heisenberg dynam- 
ics, see e.g. [38]. Improved estimates allowed for this result to be generalized, e.g., 
the case of polynomially decaying interactions was covered in [6] and the existence 
of the dynamics for the general systems considered here was proven in [27]. For 
perturbations of the harmonic system, the existence of the thermodynamic limit 
has been proven with two distinct methods, see [21] and [27]. 

An area law for gapped one-dimensional systems was proven by Hastings in [14]. 
In general, the area law conjecture states that the von Neumann entropy of the 
restriction of gapped ground states to a finite volume of size A grows no faster 
than a quantity proportional to the surface area of A. Certain aspects of Hastings' 
argument generalize to the multi-dimensional setting, e.g. a factorization property 
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of gapped ground states was proven in [19], but a proof of the area law for general 
gapped systems in arbitrary dimension remains an important open question. Some 
progress on a class of unfrustrated spin Hamiltonians appears in [22]. A review of 
these topics is contained in [29], see also [31]. 

Quantization of the Hall conductance for a general class of interacting fermions 
was recently proven in [20]. This intriguing result makes crucial use of improved 
Lieb-Robinson bounds and the methods associated with Hastings' quasi-adiabtic 
evolution. A detailed analysis of this technique, with specific regards to its implica- 
tions for perturbation theory, is the main topic of [28] . 

Finally, stability of topological order was addressed in [23, 24]. There the authors 
consider a class of Hamiltonians that are the sum of commuting short-range terms, 
such as the toric code model developed by Kitaev in [39], and proved that the 
topological order of the ground states is stable under arbitrary, small short-range 
perturbations. 

Developing a better understanding of quantum dynamics and its perturbation 
theory will be crucial in providing new insight into complex physical phenomena. 
As indicated by the number of recent generalizations and applications, the analysis 
of Lieb-Robinson bounds is a thriving area of active research which attempts to 
address this very issue. 
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